home *** CD-ROM | disk | FTP | other *** search
/ Aminet 15 / Aminet 15 - Nov 1996.iso / Aminet / misc / math / Laplace.lha / Laplace / Data / Projects.lha / physic_PI_34 < prev    next >
Encoding:
Text File  |  1996-06-28  |  1.1 KB  |  48 lines
  1. #LAPLACE
  2. #USEFLOAT=FALSE
  3. #SIMPLIFY=TRUE
  4. #ERRBND=15
  5. #EXPDSP=5
  6. #PREC=8
  7. #F1= + abs(diff(c_2,v_2)*dv)
  8. #F2=g(x)=
  9. #F3=h(x)=f(x)/g(x)
  10. #F4=matrix(
  11. #F5=vector(
  12. #F6=diff(f)
  13. #F7=diff(g)
  14. #F8=diff(h)
  15. #F9=taylor(f,0)
  16. #F10=taylor(f,1)
  17. #F11=%(run type s:startup-sequence)
  18. #F12=
  19. #F13=
  20. #F14=
  21. #F15=
  22. #F16=
  23. #F17=
  24. #F18=
  25. #F19=
  26. #F20=
  27. c_0=1e-3
  28. c_1(v_0,v_1)=c_0*(v_1)/(v_0+v_1)
  29. c_2(v_0,v_1,v_2)=c_0*(v_1+v_2)/(v_0+v_1+v_2)
  30. c_3(v_0,v_1,v_2,v_3)=c_0*(v_1+v_2+v_3)/(v_0+v_1+v_2+v_3)
  31. c_4(v_0,v_1,v_2,v_3,v_4)=c_0*(v_1+v_2+v_3+v_4)/(v_0+v_1+v_2+v_3+v_4)
  32. dv=1e-4
  33. v_0=21e-3,v_1=1.4e-3,v_2=1.6e-3,v_3=4e-3,v_4=14e-3
  34. dc_1=abs(diff(c_1,v_0)*dv) + abs(diff(c_1,v_1)*dv)
  35. eval(dc_1)
  36. dc_2=abs(diff(c_2,v_0)*dv) + abs(diff(c_2,v_1)*dv)  + abs(diff(c_2,v_2)*dv)
  37. eval(dc_2)
  38. dc_3=abs(diff(c_3,v_0)*dv) + abs(diff(c_3,v_1)*dv)  + abs(diff(c_3,v_2)*dv) + abs(diff(c_3,v_3)*dv)
  39. eval(dc_3)
  40.  + abs(diff(c_4,v_1)*dv)  + abs(diff(c_4,v_2)*dv) + abs(diff(c_4,v_3)*dv) + abs(diff(c_4,v_4)*dv)
  41. a=abs(diff(c_4,v_0)*dv),a
  42. b=abs(diff(c_4,v_1)*dv),b
  43. c=abs(diff(c_4,v_2)*dv),c
  44. d=abs(diff(c_4,v_3)*dv),d
  45. e=abs(diff(c_4,v_4)*dv),e
  46. a+b+c+d+e
  47.  
  48.